Thermophoresis of strongly charged colloidal particles 
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We present a microscopic theory for thermally driven transport in colloidal suspension of charged 
particles. For small valency, we recover known results from Debye-Hiickel approximation, which 
are in quantitative agreement with recent experiments on polystyrene particles and sodium dodecyl 
sulphate micelles. Significant modifications occur for strong charges, where the thermodiffusion coef- 
ficient is independent of the Debye length and the permittivity of the electrolyte, but is proportional 
to the valency of the macroion. 
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I. INTRODUCTION 

A thermally driven flow, or Ludwig-Soret effect, is 
observed when applying a thermal gradient to a com- 
plex fluid P, H, y, 0]. The corresponding mass trans- 
port is relevant for natural and technologicalprocesses, 
such as the global circulation of sea water [SI and the 
phase behavior of eutectic systems at solidification Q. 
In recent years, detailed experimental studies on macro- 
molecular solutions and colloidal suspensions have re- 
vealed unambiguous and often surprising dependencies 
of the Soret effect on system parameters such as salin- 
ity, surface coating, solute concentration, and molecular 
weight 0, S i, [ia, [il 111, [il, III 111, [lE [13, lill. Al- 
though the analogy to electrophoresis indicates the rele- 
vance of surface or Marangoni forces and suggests a hy- 
drodynamic treatment 1^, l^l , the physcial mechanisms 



that drive thermophoresis in liquids are poorly under- 
stood and differ from those in gaseous phases [il!, Ull . 

Diffusion theory accounts for the particle current due 
to density gradient in terms of Fick's law J = — DVn, 
where n is the number of suspended particles per unit 
volume. For a closed system the stationary state is char- 
acterized by J = and, in the absence of external forces, 
corresponds to a homogeneous density. Thermophoresis 
describes an additional flow that is driven by a tempera- 
ture gradient, resulting in the total particle current 



J = -DVtl-uDtVT 



(1) 



with the thermal diffusion coefficient Dt- Experimen- 
tally, Dt is determined by applying a temperature gra- 
dient to a uniform suspension (Vrt = 0) and by recording 
the initial current J — —uDt^T, or by measuring the 
density modulation Sn — —u^Dt / D)6T induced by a 
temperature inhomogeneity ST in the steady state J = 
[3|. The latter method gives the Soret coefficient 



St = Dt/D 



(2) 



Eq. ([T]) provides a macroscopic description for the par- 
ticle current in terms of the "thermal forces" In 
a microscopic picture, one wishes to relate the kinetic 
coefficients D and Dt to the properties of solute and sol- 
vent. For this purpose we rewrite the particle current as 



the sum of two terms }24 



J = nu — /iVn 



(3) 



The phoretic velocity u arises from the interactions at the 
solute-solvent interface and thus is a single-particle effect, 
whereas the remainder is due to the mutual forces exerted 
by the suspended particles in terms of the gradient of 
the osmotic pressure 11. The bare mobility fi — l/lGnarj) 
depends on the solvent viscosity rj and on the particle size 
a. Eq. ^ describes a non-uniform system and cannot be 
obtained from equilibrium statistical mechanics [H, |2^ . 
The stationary state J = provides the condition of 
mechanical equilibrium, i.e., when all forces acting on a 
given particle cancel. 

As an illustration, let us briefly discuss the case of an 
ideal gas. Because of the absence of interactions, one has 
u = and the osmotic pressure reads 11 = nfcsT. The re- 
sulting particle current is then J = —fikslT'^n + n^T). 
Comparing to Eq. ^ gives the Einstein relation D — 
fj-ksT as well as Dt — fJ-ks, with the Soret coefficient 
St — l/T. This simply means that, at constant pres- 
sure, the stationary density is inversely proportional to 
the non-uniform temperature and that the particles ac- 
cumulate in regions of lower temperature. Yet most col- 
loidal suspensions show a considerably stronger, positive 
or negative, Soret effect, i.e., the thermally driven current 
by far exceeds that of an ideal gas and may be directed 
toward colder or warmer regions. These deviations ex- 
press the failure of the ideal-gas picture for the solute 
and emphasize the importance of particle-solvent inter- 
actions. 

In the present paper we discuss the Soret effect in 
charged colloidal suspensions; we complete the theoret- 
ical framework given in our previous work (2^ . l25j and 
extend the electrostatic part to strongly charged parti- 
cles. In SeclUwe derive the single-particle velocity from 
the Stokes equation for the surrounding fluid, with slip 
boundary condition that is related to the anisotropy of 
the Maxwell tensor in the diffuse double layer. Sec. IIIII 
gives the thermodiffusion coefficient Dt in terms of the 
surface stress, or Marangoni force, and the osmotic pres- 
sure. In Sec. IIVI we obtain the electrostatic properties 
of the boundary layer and, in particular, the anisotropy 
parameter. Sec. |V] compare our predictions with exper- 
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imental findings and previous work, and we summarize 
our main results in Sec. IVII 



II. HYDRODYNAMICS 

The thermal forces acting on a colloidal particle and 
the resulting velocity cannot be obtained from a mechani- 
cal or thermodynamic potential. Thermophoresis is a hy- 
drodynamic problem, and the particle velocity has to be 
derived from a fluid- mechanical treatment ^2Q] . The ther- 
mally driven motion of micron or nanometer sized parti- 
cles in a viscous liquid involves small Reynolds numbers, 
i.e., inertia effects are negligible. Then the stationary 
velocity is given by Stokes' equation [H, 113, Hsl 

TyV^v = VP - V • r , (4) 

where rj is the solvent viscosity. An incompressible fluid 
satisfies moreover V • v = and stick boundary condi- 
tions at the particle surface, v = u. Besides the pressure 
gradient VP, the fluid is subject to the Maxwell stress 
tensor Tij which, contrary to the hydrostatic pressure P, 
is anisotropic and thus takes different values parallel and 
perpendicular to the particle surface. 

The force field induced by surface charges, V • T, is 
finite only within a boundary layer thickness A, which is 
given by the Debye screening length A — k^^. Through- 
out this paper we suppose that the thickness of the 
boundary layer A is much smaller than the particle radius 
a, 

A < a . 

The characteristic length scales of the normal and par- 
allel derivatives in Eq. ^ are given by A and a, respec- 
tively. Thus the forces vary rapidly in the normal di- 
rection, and much more slowly along the interface. This 
separation of length scales permits us to calculate the 
particle velocity in two steps, as schematically drawn in 
Fig. [TJ First, resorting to a ID approximation that is 
valid at distances much shorter than the particle size, we 
discuss the boundary layer and derive the surface force 
induced by the temperature gradient. In a second step 
we obtain the fluid velocity fleld at distances well beyond 
A, where the Maxwell tensor vanishes. 



A. Boundary layer 

Close to particle, the surface may be considered as flat, 
with normal and parallel coordinates z and x. Since the 
electric field E of a uniformly charged spherical particle 
is normal to the surface, the Maxwell tensor is diagonal 
with elements Tu where i — x,y, z. Then Stokes' equa- 
tion reads 
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FIG. 1: (Color online) Illustration of the boundary layer ap- 
proximation developed hereafter. Very close to the particle 
surface, a ID approximation is valid. This solution then acts 
as a boundary condition for the problem at distances well 
beyond the layer thickness A. 

with dx = d/dx, etc. The normal velocity vanishes close 
to the interface, Vz = 0; thus dz{P — T^z) = so that 
P — Tzz = Po is a constant. Inserting this relation in the 
equation for Vx and noting that, because of the thermal 
gradient, is slowly varying along the surface, we find 

vd^zVx^dx{%z-%x) ■ (5) 

The derivative on the right-hand side gives the lateral 
force per unit volume exerted on the fluid. 

Since T is finite within the boundary layer only, the 
force density dx{Tzz — Txx) vanishes beyond a distance A 
from the interface. Integrating over the normal coordi- 
nate gives the total shear stress on the fiuid accumulated 
through the boundary layer, 

-J^ dz{Tzz-Txx)^~^. 

The right-hand side gives the lateral force df acting on 
the surface element dS of the particle. 

B. Surface force on a spherical particle 

Well beyond the boundary layer, the surface can no 
longer be taken as flat and the above approximation 
ceases to be valid. At distances comparable to the par- 
ticle radius a, one has to deal with the 3D Stokes equa- 
tion, yet with modified boundary conditions. The elec- 
tric stress is absorbed in the boundary layer, and Eq. ^ 
simplifies to r/V^v = VP. It turns out convenient to de- 
fine the anisotropy parameter characterizing the bound- 
ary layer 

/"OO 

7=/ dr(r^-^i), (6) 

J a 
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C. Single-particle velocity 




FIG. 2; (Color online) Representation of the particle with 
the normal and parallel directions n and t. The polar angle 
is noted 9. 



where 71 and 7|| are the normal and parallel components 
of the Maxwell tensor, and where we have used A ^ a. 
Then the oriented surface force reads 



df „ 



(7) 



its direction is parallel to the particle surface and oppo- 
site to the 2D gradient V||7. In the case of a non- uniform 
temperature one has 



V||7 = 7tV||T 



(8) 



where we have assumed that spatial variations arises from 
the thermal gradient only and that the elements of the 
Maxwell tensor depend on temperature. 

Since heat propagation is much faster than particle mi- 
gration, the temperature field may be taken as stationary. 
Starting from the overall thermal gradient of the fluid 

VT = T,e, 

with constant T^^ the heat conduction equation for a 
spherical particle is readily solved, and the gradient at 
the particle surface reads [2^| 

V||r = $(t • Vr)t = -^T^sinM , 

where the polar angle is given by sin^^ = —e^ ■ t — see 
Fig.[2]for a definition of the unitary vectors n and t. The 
parameter 



3ks 



2ks + up 



is determined by the heat conductivities of solvent and 
particle. 



In a mesoscopic description, we may put A — > and 
consider V||7 as a surface force acting at the particle- 
fluid interface. Thus Eq. ^ reduces to the force-free 
Stokes equation 

TyV^v = VP , (9) 

albeit with modified boundary conditions that account 
for the surface force accumulated through the boundary 
layer. The stress tensor 

comprises the dissipative term or viscous force de nsity 
(T-j = "H (diVj + djVi), and the hydrostatic pressure P [2a |. 

The stick boundary conditions of hydrodynamics, im- 
posing a continuous velocity at the interface, vj^,^^ = u, 
cease to be valid in the presence of surface forces. The 
fluid being incompressible, the normal component is con- 
tinuous, 

n-v|^^^ = n-u, (10) 

whereas the parallel velocity in general exhibits a jump. 
A second condition is obtained by noting that there is 
no net external force acting on the system consisting of 
the particle and the boundary layer. Thus the integrated 
normal force outside the boundary layer vanishes [20|, 



dScr • n = 



(11) 



The third condition involves the surface force ([5]), which 
plays the role of the source field and imposes a finite shear 
stress (H, 



t • (cr • n-h V[|7) = 



(12) 



It turns out that this condition results in a finite slip 
velocity. 

The source term (IT2]) differs from that of electrophore- 
sis, where the Stokes equation in the boundary layer 
reads rjd^Vx — eEcxtd'^ip, with the external electric field 
Ecxt and the local potential ip. Integrating this equation 
twice, with stick boundary conditions, gives the meso- 
scopic slip velocity vs = — eiJoxtV'o/^ as a function of the 
surface potential V'o, resulting in t- (v — u)|r=a = vs 
In the case of thermophoresis, integration of Eq. (O re- 
sults in the surface stress V||7; the slip velocity is ob- 
tained only after solving Eqs. (P))- p^ . 

In the laboratory frame the particle moves at a velocity 
u = ue^. Transforming to the reference frame in which 
the particle is at rest, we have u = and v(r) = v(r) — 
u. The solution of Stokes' equation at small Reynolds 
numbers in spherical coordinates v = Vrii+igt reads psj 



i;^ = -ucose* ( 1 - 2a- + 2I3-. 

r r 



ve ~ u sm f 



a a" 
1 — a p — 



(13a) 
(13b) 
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where 9 is the polar angle with respect to the x axis 
and the radial and tangential unit vectors n = r/r and 
t = dn/dO satisfy ea, — cosOn — smQt. This flow field is 
related to a non-uniform hydrodynamic pressure 



P(r) = Po + a 



2r]ua 



■ cosf 



D. Mobility 

Eq. ([7]) accounts for the balance of the stress on the 
fluid and the force di acting on a surface element dS of 
the particle. Inserting the expression of V||T and sepa- 
rating the dependence on the angle 9, we have 



The parameters u, a, (3 are determined from the solulu- 
tion of Stokes' with the boundary conditions (fT0| - (fT2)) . 

The first two of these conditions involve the fluid veloc- 
ity and stress only. In the particle-fixed frame the normal 
velocity vanishes, Vr\r=a — 0, resulting in 1 — 2a + 2/3 = 0. 
The total stress at the interface can be written as cr • n = 
narr + 'tcTre, with the entries of the dissipative part in 
shperical coordinates [28l | 



27? — 

dr ' 



dv0 
dr 



r 



Integrating pT|) over a sphere just outside the boundary 
layer gives the relation 1 — 5a + 2(3 = 0. One readily 
obtains the amplitudes of the velocity field varying with 
distance as 1/r and l/r^, respectively, 



a = 0, 



(14) 



The remaining condition (fT2|) relates the stress a'^g and 
the surface force and determines the particle velocity 



' 377 



T 

J. rp 



(15) 



Finally, taking the back transformation v(r) = v(r) + u 
yields the fluid velocity in the laboratory frame 



2Tra^jTKT^ / d{cos9)sin9t 



The integral takes the value 
force 



-^Gx, resulting in the total 



f = _ 7ra2^^^Vr . 
o 

A phoretic mobility may be defined as the ratio of the 
one-particle velocity and the force, u = /Xphf . Comparing 
with pT|) gives 



1 



(19) 



This expression differs from the mobility fi that de- 
scribes a particle subject to an external force Foxt or the 
Brownian force leading to diffusive motion. The solution 
of the Stokes equation with such an external force satis- 
fies stick boundary conditions v|^^^ — u and the force 
balance Foxt + / dScr • n = 0, resulting in 



3 R 1 



Up = 



the excess hydrodynamic pressure P(r) — Pq ~ 1/f^, and 
a velocity field ^/'p that decays with distance as ~ l/r. 



v(r) 



sin 9t 



cos Wn 



(16) 



The slip velocity, i.e., the velocity change through the 
boundary layer, is given by the tangential component 
veia), 



V5 — -usin9t. 



Its maximum value ^Ti at 



= ^ exceeds the particle 
velocity, which means that the particle and the fluid move 
in opposite directions. 

Eq. provides the relation the phoretic velocity u 
and the applied temperature gradient VT = TxQx, 



u = -CVT , 
with the transport coefficient 



C = 



d'y 
3^dT 



(17) 



(18) 



This result displays how the phoretic velocity arises from 
the competition of particle-solvent interactions that drive 
the particle, and the fluid viscosity. 



III. KINETIC COEFFICIENTS 

The results of the previous section are now applied to 
derive general expressions for the transport coefficient D 
and Dt- We consider a suspension of interacting parti- 
cles in an inhomogeneous temperature field r(r), where 
both the density n and the osmotic pressure 11 are, in 
general, spatially varying functions. Such a system can- 
not be dealt with in terms of equilibrium statistical me- 
chanics; its steady state and the relaxational kinetics are 
described in terms of the particle current Eq. After 
expressing both terms by microscopic quantities propor- 
tional to Vn and VT, we compare to Eq. ([1]). 



A. Osmotic pressure 

The particle current contains, besides the single- 
particle velocity u, a term proportional to the pressure 
gradient VII. The latter depends on both density and 
temperature modulations. 



vn = nTVr + n„Vn 



(20) 
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with the shorthand notation IIt = dH/dT, etc. 

For a non-uniform system with a pair interaction V ^ 
the pressure is calculated by explicitly writing the forces 
on a particle at position r; proceeding as with the Clau- 
sius virial function one finds the sum of the idcal-gas pres- 
sure nksT and the two-particle forces iij = —'ViV{ri — 

r,) m, 

n(r) = n(r)fcBr(r) + ^ ' ^^^^^"^ " ^»)^ ' 

Assuming the pair interaction V{rij) to be isotropic 
and sufficiently short-ranged, one has with the density- 
density correlation function g{r,r') 

n = nksT-^ J dv'giryy^^ . 

This relation is a useful starting point for a perturbative 
evaluation of the diffusion current in terms of tempera- 
ture and density gradients. 

So far we have avoided the use of equilibrium ther- 
modynamics; the above relations involve mechanical and 
hydrodynamic forces only. Since D and Dt are linear- 
response coefficients, they may be evaluated in terms of 
equilibrium statistical mechanics, and (...) may be con- 
sidered as the canonical average with the Gibbs distri- 
bution. This means in particular that the density and 
temperature fields are replaced by their mean values n 
and T. 

For dilute and weakly interacting systems, a widely 
used approximation consists in expanding the pressure 
in powers of the particle density, 

ksT {n + Bn^ + . ..) , (21) 

and truncating at quadratic order, with a coefficient 



B 



1 - e 



-V{r)/kBT 



(22) 



that still depends on temperature (29[ . Most relevant pair 
interactions are repulsive, such as the screened electro- 
static potential in an electrolyte and the entropic short- 
range repulsion due to polymers grafted on colloidal par- 
ticles. Then the virial coefficient is positive, B > Q and 
the interactions enhance the pressure. 



B. Thermal diffusion coefficient 

Inserting u and VII in Eq. ([3]) and comparing to 
Eq. H]) yields 



D = ^n„ , and Dt = -Ht + C 



(23) 



Together with the virial expansion Eq. (HI]), the transla- 
tional diffusion coefficient takes the form 



which does not depend on the temperature gradient and 
is insensitive to surface forces. A different picture arises 
for the remaining coefficient. 



Dt = ^lkB {l + Bn + hTBt 



C 



(25) 



which depends explicitly on the temperature derivative 
Bt and on the surface force parameter C. 

For an ideal gas both the particle-solvent potential and 
the pair interaction vanish, resulting in C = and B — 0. 
One readily recovers the transport and Soret coefficients 
D — nksT, Dt = jJ^ks-, and St = 1/T. 



IV. THE ANISOTROPY PARAMETER 

The above expression for the transport coefficient de- 
pends on the charged double layer through the parame- 
ter 7 which characterizes the anisotropy of the Maxwell 
stress tensor. The electrostatic properties of the bound- 
ary layer are closely related to the distribution of the 
mobile co-ions and counter-ions of the electrolyte. In 
three dimensions, the resulting Poisson-Boltzmann (PB) 
equation cannot be solved analytically. Approximate so- 
lutions exist for two particular situations. First, for weak 
charges, the potential e5'(r) is small as compared to the 
thermal energy and the PB equation can be linearized. 
Second, for particles that are significantly larger than 
the Debye length, the potential can be evaluated as an 
asymptotic series in terms of the small parameter 

(26) 

Since a small radius a < n^^ implies rather a weak 
charge, and large valencies are possible for sufficiently 
large particles only, the two approximations cover in fact 
most physical systems. In both cases, the resulting po- 
tential ^E* is an effective quantity that contains, besides 
the potential energy, the entropy of the mobile ions. 

The limiting cases of weak and strong charges are ex- 
pressed by the ratio of the bare valency Z and its effective 
value 



Z* 



4(1 + Ka)a 



(27) 



that is given in terms of the particle size a, the Debye 
screening parameter 



(28) 



and the Bjerrum length £b = e^/(47refcBT), with e the 
dielectric constant of the solvent (e/eo — 78 for water). 



A. The Poisson-Boltzmann equation 

In the mean-field approximation, the electrostatic po- 
tential of a charged particle satisfies the PB equation 



D = ^ikBT{l + 2Bn+ ...) , 



(24) 



^2 X / N 2noe . , e4'(r) 

V^^'(r) = — — sinh — 

^ ' e kuT 



(29) 
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This equation has to be solved together with the follow- 
ing boundary conditions. First, the potential vanishes 
at infinity, linir^oo '^{r) = 0, and second, the electric 
field at the surface of the colloid is related to the surface 
charge a according to 



dr 



a 

e 



(30) 



However, since the PB equation is non-linear, one can 
only find an approximate solution in spherical geome- 
try. We discuss below the two relevant limits of weakly 
charged and highly charged colloids. To proceed, it turns 
out convenient to define the reduced potential 



V{r) 



that obeys the differential equation 
V^y — sinhy 



(31) 



(32) 



B. Debye-Hiickel approximation 

We first recall the electrostatic properties for weakly 
charged systems. In this case, the reduced potential is 
small everywhere, y <C 1. Then Eq. ([5^ can be lin- 
earized, V^yDH(r) — K^yDH(r), and solved by the well- 
known Debye-Hiickel potential 



2/DH(r) = 



1 + Kfl r 



(33) 



This solution ceases to be valid if the surface potential 
2/DH(a) approaches unity, i.e., if the valency Z is of the 
order of the effective value Z* . 



and vanishes for n> 1 



dr 



(36) 



The multiple scale approach [30| allows to systematically 
match the known ID solution close to the particle surface 
to the screened potential at distances well beyond the 
Debye-Hiickel length. This is achieved by introducing 



Xi 



I /r 
e 



X2 



(37) 



where the "fast" variable xi accurately describes the 
problem in the vicinity of the surface, in addition to 
the "slow" variable X2 which accounts for the spheri- 
cal symmetry of the problem far from the particle. The 
derivative thus becomes d/dr = a^^{e^^di + 82), with 
di = d/dxi, and the Laplace operator reads in terms of 
this new set of variables 



V2 = df + 2e( 8^82 + -92) dl + -82 

\ Xi J \ X2 

Finally we expand the S — sinh y according to 
5=.5o + (et/i+e^ya + .-OCo 



(38) 



I 



■5o 



with 5o = sinhyo and Co = coshj/o- Inserting these rela- 
tions in the PB equation ([5^ , we obtain a set of coupled 
equation for ?;„ that can be solved by identifying terms 
of the same power of e. We give explicitly the equation 
at order e*^. 



dlyo = sinh 7/0 



(39) 



and that at order e^. 



C. The multiple scale method 

For highly charged particles the nonlinearities in the 
PB equation are essential and result in a much more in- 
tricate problem. Still, an asymptotic expansion 



yir) = yoir) + eyi(r) + e'^y2ir) + 



(34) 



can be derived explicitly if the parameter e is small, i.e., 
if the Debye length is smaller than the particle radius. 
Note that the coefficients 2/„ in still depend on e, 
i.e., we have not yet defined how to construct the series. 

The usual boundary conditions of electrostatics impose 
that the potential vanishes at infinity, limr^oo Vnif) = 0, 
and satisfies (|30p at the particle surface. Since Eq. (|30p is 
independent of Ka, the derivative involves the term ?i = 
only 



dr 



ksTe ' 



(35) 



8lyi + yi coshyo = — 8iyQ + 28i82yo ■ (40) 

X2 



1. Zeroth-order solution yo 

At lowest order in e the differential equation ([39]) can 
be integrated as (3l| 



diyo 



-2 sinh 



yo 



where we have used the condition yo ^ and diyo ^ 
as xi 00. This equation can be integrated once more 
and one finds [sjj 



yo 



2 In 



l + /(x2)e-"^ 
I-/(x2)e-^i 



where the unknown integration function /(X2) remains 
to be determined. 



7 



The point is now that we artificially introduced an ex- 
tra variable in the problem. The equation of order and 
the boundary conditions at r = a and r ^ oo are not suf- 
ficient to obtain the function of two variables ^0(2:17 2^2)- 
Following [32] , we thus require that the non homogeneous 
part of Eq. (j40|) for the first order solution decays faster 
than exp[— a;i]. Explicitly, this requirement reads 

X2d2f{x2)+f{x2)=0 

so that /(X2) — go/x2, with 170 constant. 

Coming back to the original variable, the lowest order 
term of the expansion reads 



Voir) = 21n 



1 + gjr) 

1 - 9ir) 



with 



air) = go-exp[-K(r- a)] 



(41) 



(42) 



The integration constant go is finally obtained from the 
boundary condition (|35p 



90 



+ A 1 



2Z 



(43) 



To leading order in 1/ (na) , the parameter p may be writ- 
ten as p = /£gc, i-e., the ratio of Debye and Gouy- 
Chapman lengths k^^ and iac = l/{2TTa£B)- 



4- Higher order corrections 

Higher order terms in the expansion can be evaluated 
following the same method. We give below only the first 
order correction, which is obtained by solving Eq. (|40p to- 
gether with the requirement that the non homogeneous 
part of the equation for y2 does not contains terms pro- 
portional to exp[— xi]. We obtain 



yi(r) 



2g(r) 



1-giry 



-g{rf + k 



The constant k is obtained by enforcing the boundary 
condition at the surface of the particle (dyi/dr)^^^ — 0. 
Explicitly, we find 

(l + 2e).g^,(3 + 4e).gg 

(l + .9o')(l + e) 

Expanding this expression at lowest order in e, one re- 
covers the result of Ref . [32] . Incidentally, we note that 
the effective charge obtained from the behavior of the 
potential in the limit p 3> 1 is now corrected to become 
Z* = jf-(| + na), as pointed out by Trizac and collabo- 
rators [33j ]. This correction is, however, irrelevant in the 
limit na ^ 1 considered here. Contributions of quadratic 
and higher order in e can be obtained by repeating the 
same procedure. 



D. The anisotropy parameter 



2. Small valency Z Z* 

If p <C 1, one has go — p/2 and the solution ([iT 
reduces to 



ZIb e-'^^'-") 
1 + Ka r 



(44) 



As expected, we recover the Debye-Hiickel solution for 
weakly charged colloids. 



In view of Eq. ^ , we now evaluate the Maxwell stress 
tensor which is given by the components of the electric 
field E = -V*, 



1 



Since the electric field of a spherical particle is normal to 
the surface, the off-diagonal or shear components vanish; 
the diagonal elements are finite, the normal one being 
positive and the two parallel to the surface negative. 



3. High valencies Z Z* 

In the limit p 1, the solution y^ exhibits quite 
different behaviors at short and long distances. In the 
latter case, one recovers well beyond the Debye length, 
K.(r — a) 3> 1, the screened potential 



ZHb e 



-K{r — a) 



1 



(45) 



In the limit Ka 3> 1 the particle surface may be con- 
sidered as flat, and the anisotropy parameter Q reads 



7 = e / drEirf . 

J a 



Inserting j/o one obtains 7 = e{kBT/e)'^ J dr{dyo/dr)^. 
Replacing the valencies by the charge number densities 



albeit with the effective valency Z* . Close to the surface 
of the particle, K(r — a) <C 1, the potential simplifies to 
a logarithmic law 



yo{r) — — 21n/t(r — a) . 



(46) 



47ra2' 



Ana^ nil 



we find 



(47) 
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In the two Hmiting cases with respect to the ratio a ja* 
one has 



7 



_ ( 2kBTa^/a* for o" < cr* 



2kBTa 



for a* a 



(48) 



The effective charge density a* depends on temperature, 
whereas the bare value a is constant. In an inhomoge- 
neous system, T(r) is slowly varying along the interface, 
and a* and 7 are well defined local quantities. 

The derivative 7t is calculated by rewriting the effec- 
tive charge density as a* — {k/t:£b) and by noting the 
dependencies of Debye and Bjerrum lengths, fi~'^ cx V eT 
and cx l/(eT), and of the dielectric constant e ^34i] 

dine 
" " ~dh[f ■ 

Thus one finds the relation da* /dT = (1 - a)a* /2T, and 
finally obtains 

At room temperature the parameter a takes the value 
a = 1.4. 

In the limit of weak charges one has 77- — /T){1 + 
a) and thus recovers the expression given previously 
in M 



7T = fci3(l + a)— for cr< cr* , (50) 
cr 

whereas in the opposite case we find 

7t = 2/cBcr for a* -^a . (51) 

These expressions for the parameter 77^ constitute a main 
result of the present paper; with Eqs. and they 
provides the explicit form of the thermodiffusion coeffi- 
cient Dt- 

V. DISCUSSION 



A. Dominant driving force 

It is instructive to estimate the relative magnitude of 
the two contributions to Dt, that is the ratio 



C 



47ra 7t 



(52) 



This quantity exceeds unity for most experimentally rel- 
evant systems. Thus for micron-sized polystyrene beads 
with typical parameters a ~ \ ^m, cr « 10~^ e/nm^, 
~ 10 nm, and the viscosity and permittivity of wa- 
ter [1^, one finds the value C/fiks ~ 10'^. Even for 
solutes in the nanometer range, this ratio may be large; 
for instance in the case of sodium dodecyl sulphate (SDS) 
micelles of valency Z « 30 and radius a « 3 nm 

[lal, the 

resulting ratio C / fjiks lies between 1 and 10. Note that 
we have used here ~ 1, which holds for most colloidal 
suspensions. 

Similar arguments hold for polymer solutions, where 
the monomer-solvent interaction is at least of the order 
/cbT. The phoretic coefficient C ~ kB/rjo. scales with 
the inverse monomer size a, whereas the mobility /i ~ 
\/{rjR) involves the gyration radius R ^ N'^a, resulting 
in C/fikB ~ TV > 1 [2I. 

Thus we are led to the important conclusion that for 
solutes larger than a few nanometers, thermal diffusion 
is in general driven by surface forces, i.e.. 



Dt^C 



(53) 



This statement is confirmed by the experimental obser- 
vation that hardly any system satisfies the ratio Dt /D = 
1/T expected for diffusive thermal transport. In physical 
terms this means that the ideal-gas or entropic contribu- 
tion /ifcs is much smaller than the surface force term C 
that arises from solute-solvent interactions. Note also 
that 7t is always positive for the electric double layer, 
which means C > and a thcrmophoretic current to- 
ward colder regions. 



The preceding sections provide an explicit expression 
for the thermally driven particle current 

Jth = -uDtVT , 

which according to Eq. consists of two contributions; 
the first one, n^kB'^T, describes a diffusive current due 
to the non-uniform stationary state in a temperature gra- 
dient, whereas the second one, nCVT, arises from sur- 
face forces. The latter contribution is proportional to the 
parameter 77^ calculated above. 

In this section we discuss the relative magnitude of 
these contributions, their dependencies on physical pa- 
rameters of the solute and the solvent, and compare our 
results with experimental findings and previous theoret- 
ical work. We mostly consider the dilute limit, where 

Dt = ^J.kB + C . 

Interaction effects that occur at higher densities are ad- 
dressed separately. 



B. Particle size 

In the range of parameters considered in this paper, 
Ka ^ 1, the parameter cr* is independent of the particle 
size, and so are the surface energy 7 and its derivative 
7t. As a consequence, the phoretic coefficient C varies 
linearly with the particle size, 

C oca . (54) 

The stationary velocity u = —CVT is thus proportional 
to a, whereas the diffusive contribution u = —^kB^T 
varies as 1/a. 

C. Surface charge and Debye length 

For the case of weakly charged surfaces a <^ a* , one 
finds from Eq. that the phoretic coefficient C varies 
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FIG. 3: Reduced phoretic coefficient {3rjC) / {2kBO-* ^a) as a 
function of tlie reduced ctiarge density a/a* . The beliavior of 
tlie pfioretic coefficient switclies from quadratic at iow cfiarge 
density to iinear at liigh cfiarge density. 



quadratically with the charge density 



C = 



12?]Te 



{l + a)(7^aK ^ {(7^(7*) (55) 



Moreover, C is proportional to the particle size a and to 
the Debye screening length k^^. Both dependencies on a 
and X = are confirmed by experiment, as discussed 
below for the Soret coefficient. 

In the opposite limit of strongly charged particles, the 
phoretic coefficient 



C 



2^ 
3?7 



■era ((T 3> <y*) 



(56) 



varies linearly with the surface charge density and with 
the particle size, but is independent of the Debye length 
and the dielectric constant of the electrolyte. The depen- 
dence with the reduced charge density a /a* is shown on 
Fig. [31 



D. Soret coefficient 

For the case of small valencies cr <C cr*, the Soret coef- 
ficient reads 



1 



T 



2 efcsT 



(57) 



As discussed below, the ideal-gas value 1/T is small for 
particles larger than nanometers and may be safely ne- 
glected. The Soret coefficient is given by the remaining 
charge term. 

The power laws St oc and St cx; A, as first predicted 
in [25j with Z = Air a? a, have been confirmed recently in 
experiments on polystyrene beads in the parameter range 
a = 100 ... 550 nm and A = 2 ... 13 nm [H, HI]. These 
dependencies may be combined to an interaction volume 
47ra^A, i.e., a spherical shell of thickness A and mean 



radius a, thus giving clear evidence for surface forces op- 
erating in a boundary layer at the particle surface. 
In the limit ct 3> cr* , the Soret coefficient reduces to 



1 -f $47rcra2 1 + 



T 



T 



(58) 



For large valencies one has the particularly simple form 
St = iZ/T. 



E. Interaction effects 

Now we turn to higher particle densities where the first 
correction of the virial expansion is no longer small. Al- 
though the surface-force driven coefficient C does not 
depend on density, the diffusion coefficient varies as 
D — iiksT^l + 2nB). Thus for a repulsive pair inter- 
action the Soret coefficient decreases as a function of the 
density 



STin) 



ST{n = 0) 
1 + 2nB 



(59) 



As shown in Appendix \^ the virial coefficient B de- 
pends on the Debye length A, which in turn depends on 

— 1/2 

the electrolyte strength as A oc Ug . At high salin- 
ity, one has a » A, resulting in strong screening of the 
electrostatic interaction and the excluded volume of hard 
spheres, i.e., 2B = ^{2a)^. In the dilute limit or at high 
salinity, the Soret coefficient is proportional to A and thus 
varies as 



.-1/2 



(low density) 



As salinity decreases, the Debye length increases, and so 
does B. If the screening length becomes comparable to 
the particle size, a ~ A, the virial coefficient varies as 
B ~ rig ^^'^ . As a consequence, for a dense suspension in 
a weak electrolyte, the above behavior is modified by the 
virial coefficient, and the Soret coefficient increases with 
salinity, 

St ~ no (high density) . 

Both dependencies on n and no have been, at least 
qualitatively, observed for SDS micelles in aqueous so- 
lution [l3]. The inverse coefRcient 1/St coefficient has 
been shown to vary linearly with density, and its slope to 
decrease as salinity augments, like the factor 2B above. 
Moreover, Fig. 1 of Ref. [T3| reveals that in the dilute 
limit, increasing the salinity reduces the Soret effect, 
whereas for the dense micelle suspension St increases 
when adding salt. These experimental findings agree al- 
most quantitatively with our theoretical expressions. A 
more detailed comparison is not an easy undertaking, 
since both the size and the charge of the micelles slightly 
vary with the electrolyte strength. Also, part of the data 
are in the range a X whereas the present work relies 
on a 3> A. 
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F. Negative Soret effect 

The charged double layer is predicted to result in a 
"normal" Soret effect St > 0, that is, the particles 
drift in the direction opposite to the thermal gradient 
and thus accumulate in colder regions of the sample. 
There are, however, systems that show a more complex 
behavior. Aqueous suspensions of maghemite particles 
present an inverse Soret effect St < 0, whereas the nor- 
mal behavior has been observed in organic solvents ^1^- 
Also, the Soret coefficient of macromolecular solutions of 
lysozyme [T3| or DNA [l^ changes sign as a function of 
temperature and salinity. Especially these latter exam- 
ples strongly suggest that the observed transport coef- 
ficient Dt comprises contributions of different physical 
origin and of opposite sign. 

Fig. 3a of Rcf. 18] displays particularly instructive 
data for the Soret coefhcient of polystyrene beads as 
a function of the Debye length A at different tempera- 
tures. St varies linearly with A, as expected from (f57| . 
yet does not vanish at A = 0; the finite intercept is nega- 
tive at 6°C, increases with temperature, and changes sign 
at about 20° C. These observations give evidence for an 
additional contribution S'rp that is insensitive to charge 
effects but strongly temperature-dependent. A similar 
behavior has been reported for macromolecular suspen- 
sions [Ulil. 

These experimental findings indicate that, besides the 
Maxwell stress tensor of the charged double layer, ad- 
ditional forces operate at the particle surface [H, \3T\ . 
For micron sized solutes, solid-fluid interface tension and 
grafted surfactant molecules contribute to the anisotropy 
parameter 7. A very rough estimate for the tension 
may be obtained from the phcnomenological law 7/ = 
7o(l — T/Tq) for its temperature dependence. With pa- 
rameters 70 ~ 10 mN/m and Tq « 10^ K, the force 
parameter dji / dT would by far exceed the double-layer 
contribution. Yet one should be aware that the actual sit- 
uation at solid-fluid interfaces, especially at polar ones, 
is significantly more complex than suggested by the sim- 
ple relation dji/dT = —jq/Tq. Moreover, because of 
the short range of dispersion forces, the resulting stress 
anisotropy occurs in a narrow layer of about one nanome- 
ter thickness. Macroscopic hydrodynamics with constant 
viscosity may well not apply at such short distances. 



G. Comparison with previous work 

Finally, we want to compare our analysis with previous 
work. In the spirit of Refs. [l], our description of 
thermophorcsis relies on hydrodynamics at low Reynolds 
numbers. In this approach, the temperature dependence 
of the Maxwell tensor results in a shear stress parallel to 
the surface, and the force acting on the particle is given 
by the stress accumulated through the boundary layer. 
Solving Stokes' equation with the appropriate boundary 
condition, one obtains the velocity of the particle and the 



flow of the surrounding fluid. 

A different route has been taken by several authors in 
the past few years [H, [H, [11]. In those works, the 
driving force is related to the gradient of the effective 
energy W of the charged double layer, 

F = -VW=-^VT. (60) 

This expression relies on the assumption that W may be 
treated like an external potential. The phoretic velocity 
is obtained from the steady state where the "external" 
force F and the Stokes drag F5 = —Qirrjau cancel, result- 
ing in u = —iJ,{dW/dT)'VT and in the thermodiffusion 
coefficient 

Comparison with Eqs. ([57]) and ([55)1 reveals several dis- 
crepancies which we discuss briefly. 

At small valencies the double-layer free energy and the 
integrated anisotropy parameter 7 turn out to be identi- 
cal, 

W = ina^j {Z < Z*) . 

Eq. (|B3p of Appendix |B] leads to the Soret coefficient 
St = (l + in^a^eBa^X) /T. Comparing with Eq. ^ 
shows that the charge terms obtained from hydrodynam- 
ics and from the external-force picture differ by a factor 

2 2 2ns + 

The parameter ^ accounts for the local modification of 
the temperature gradient due to the different heat con- 
ductivities of solvent and particle. For ns = up one 
has ^ = 1 and a constant gradient VT = Tx^xi whereas 
Ks <C Hp flattens the temperature profile around the 
particle and thus significantly reduces thermal transport. 
The latter case is relevant for suspensions of crystalline 
particles, where the large heat conductivity indeed sup- 
presses thermophorcsis. 

The remaining factor i occurring in the hydrodynamic 
treatment is related to the shape of the particle. For 
a prolate particle oriented along the thermal gradient, 
this factor increases, whereas it is reduced for an oblate 
shape. This dependency could be relevant for suspended 
rods, such as DNA with a persistence length exceeding 
the size of the molecule. Although it misses the numerical 
prefactor ^/2, the picture of an external force provides, 
for weak charges, qualitatively correct results. 

An additional discrepancy arises in the limit of high 
valencies. Comparing the charging energy W and the 
anisotropy parameter Ai:a?^ — 2ZkBT reveals, that they 
differ by a logarithmic factor InZ. Thus using W for 
strongly charged particles results in overestimating the 
Soret coefficient by a factor ~ InZ. 

From a more formal point of view, Eq. (j60p suffers from 
two major shortcomings. First, since the entropy S" of a 
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classical system is defined with respect to an arbitrary 
reference value Sq only, the free energy F = U—T{S—So) 
comprises an arbitrary constant So, and its the absolute 
value of dF/dT — Sq — S does not carry physical infor- 
mation. This problem persists for the charging energy 
W. Since W is defined with respect to the state where 
the excess charges are at infinity; yet this reference state 
comprises the entropy of the mobile ions, which results 
in an arbitrary contribution to the derivative dW/dT. 

A second shortcoming occurs when applying (|60p to 
polymer solutions and dense colloidal suspensions. Since 
the thermally induced surface stress ([7]) is treated as an 
external force [l^, [H, the resulting flow field 
v(r) comprises, contrary to Eq. p6p . a term ^ 1/r 
that would result in long-range hydrodynamic interac- 
tions [131. As a consequence, for a particle suspension 
of volume fraction 0, one would find a reduced value 
(1 — C,4>)Dt with C ~ 1 [39I, fi^l, which is, however, an 
artefact of the external-force picture. An even more se- 
rious discrepancy arises for polymers, where ([6T|) leads 
to underestimate the thermophoretic current by a factor 
N~'^ , with the number of monomers A'^ and the expo- 
nent v ^ 5 • ■ • I [HI • The correct treatment of the sur- 
face stress gives v(r) ~ l/r^ as in ((T6)) . i.e., there are no 
long-range hydrodynamic interactions. 

Still, the Debye-Hiickel approximation (|B3[) provides 
a correct description for dilute suspensions of weakly 
charged particles [H, [H, [s^ . Given the uncertainty of 
the surface charge density, the missing numerical prefac- 
tor is of little practical consequence. 



VI. SUMMARY 

Starting from Stokes' equation with a force term aris- 
ing from the electric stress in diffuse double layer, we have 
obtained the thermodiffusion coefficient Dt for charged 
colloidal particles. We brieffy summarize our main re- 
sults. 

(i) In the limit of small valencies we confirm previous 
work obtained in Debye-Hiickel approximation. Novel 
results arise for strong charges, where the Sorct coeffi- 
cient St = S.Z/T is proportional to the ratio of thermal 
conductivities ^ and the valency Z only. 

(ii) The coefficient Dt oc a is proportional to the par- 
ticle size, for weak charges it is moreover linear in the 
Debye length A. For strong charges, however, Dt does 
not vary with A; in this range the Soret effect does not 
depend on the salinity nor on the permittivity of the elec- 
trolyte. 

(iii) Comparing the charging energy W and the inte- 
grated anisotropy parameter 47ra^7 of the Maxwell stress 
tensor, we find that these quantities coincide for weak 
valencies, which explains the qualitatively correct re- 
sults obtained from the simple external-force picture of 
Eq. (j60p . Yet significant discrepancies arise in the oppo- 
site case of high valencies, where 47ra^7 = 2ZkBT results 
in linear behavior St Z . In this range W carries a log- 



arithmic factor ~ In Z. The resulting ambiguities would 
alter the Soret coefficient obtained from (pT|) . and illus- 
trate the need of a rigorous treatment. 
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APPENDIX A: INTERACTING PARTICLES 

In view of the expressions for the diffusion coefficients 
D we evaluate the virial coefficient B for charged parti- 
cles in a simple approximation. A "particle" consists of 
one colloidal macro-ion and its diffuse cloud of counteri- 
ons and co-ions. For repulsive potentials, the coefficient 
B accounts for an excluded volume and is expressed in 
terms of an effective particle radius vq that is larger than 
a. 

At distances r well beyond the Debye length, the pair 
interaction V{r) of two colloidal particles is given by the 
electrostatic energy of the first one in the Debye-Hiickel 
potential of the second. 



Vir) = ksT 



(l + a/A)2 r 



£g-(r-2a)/A 



(Al) 



Here Z indicates the bare valency if Z ^ Z* , and Z = Z* 
in the opposite case. 

Since the virial coefficient B cannot be evaluated as 
it stands, we resort to an approximation that consists in 
replacing the integral in (|22p by an excluded volume 



27r 
B = — r, 



' 



where ro is the effective particle radius. 

Since the particle may be considered as a rigid body, tq 
is always larger than a, thus accounting for the repulsive 
interaction. We define the effective radius as the distance 
where the pair potential is equal to the thermal energy, 

V{ro) = keT . 

Indeed, the neighbors hardly penetrate in a region of size 
ro about a given particle, whereas they move rather freely 
in the remaining space. From the expression of the pair 
potential V{r) one finds 



ro = 2 (a + xA) 



(A2) 



where the quantity 



X = 111 
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closes an implicit equation for rp . In the absence of inter- 
actions, or at very high salinity (A — > 0), the excluded- 
volume radius is twice the particle size, vq — 2a. The 
enhancement due to the electrostatic repulsion is of the 
order of the Debye screening length. With typical pa- 
rameters for SDS micelles, Z — 40, A = 1 nm, a = 3 nm, 
one finds a logarithmic prefactor of the order % 3. 



APPENDIX B: THE DOUBLE-LAYER FREE 
ENERGY 

The mean-field free energy of the charged double layer 
is usually given in terms of the "charging energy" , i.e., 
the work required to assemble the charge Q ~ Aira^ea on 
the particle from infinity, 



(Bl) 



In terms of the reduced potential yo given in ((4T|) , we find 
with Ana^a = Z 



W^ksT [ dZ'yoia) . 
Jo 

Elementary integration gives 

p + p P 



W = 2Z/crT In 



P-P P 



(B2) 



with the shorthand notation p — y'jp + l — l. In the lim- 
iting cases with respect to the ratio of bare and effective 
valencies this expression simplifies to 

r 2{ZyZ*)kBT for 

^ 2ZkBT\n{4:Z/Z*) for Z* < Z ' ^ ' 
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